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We study the two-point correlator of a modified Confined-Coulomb transition order parameter in 
four dimensional compact U{1) lattice gauge theory with Wilson action. Its long distance behavior 
in the confined phase turns out to be governed by a single particle decay. The mass of this particle 
is computed and found to be in agreement with previous calculations of the 0^"'' gaugeball mass. 
\^ , Remarkably, our order parameter allows to extract a good signal to noise ratio for masses with low 

' statistics. The results we present provide a numerical check of a theorem about the structure of the 

' Hilbert space describing the confined phase of four dimensional compact U (1) lattice gauge theories. 

(N : 

^ ' INTRODUCTION 



Quantum Chromo-Dynamics (QCD), the present theory for strong interaction, is formulated in terms of quark and 
^ giuon fields. The latter are carriers of the SU{i) symmetric colour interaction. At low energies, however, the only 
existing particles are hadrons. To explain this fact, one has to suppose that quarks are confined inside hadrons. The 

■ attempt to understand this obscure property of quarks, confinement, has led to study simplified models, defined on 
^ . a lattice, but still exhibiting confinement Q . The simplest among such models are Ising systems on two and three 

' dimensional lattices. They possess global discrete symmetries. Their generalisation to local symmetries are known as 
(«•->) ■ gauge Ising models and still exhibit confinement 0, Q ■ ^ further step to approach QCD is to consider models with 
continuous gauge symmetry in four dimensions (4D). 
. Among them is compact U{1) lattice gauge theory (U(l)lgt) 0. This model has two phases: a confined phase 

■ at strong coupling and a Coulomb phase at weak coupling. The order and the location of the phase transition that 
. separate them depend on the action form and has been a long standing subject of debate [MS Sill- The idea that 
' monopoles could play a crucial role in the description of confinement appeared after the inspiring work of Polyakov 
, in three dimensions and the contemporary conjecture, known as dual superconductivity (DS), by Mandelstam 

and 't Hooft about confinement in QCD Jii, JTj. This conjecture states that confinement could be described by a 
D . mechanism similar to that responsible of superconductivity. In the confined phase, the vacuum should be filled with 
' monopoles in the same way a superconductor ground state is filled with Cooper pairs. In the case of U(l)lgt a finite 
^ ■ density of monopoles was discovered in the confined phase and vanishing in Coulomb phase [l^ . 

The introduction of a well defined order parameter for the Confined - Coulomb transition in U(l)lgt [l^ 113 . Hi! 
proved, eventually, the condensation of monopoles at the transition. Using the Wilson action, this transition can be 
shown to be of first order (with very long correlation length) and located, in the thermodynamic limit, at /3c ~ 1.011 
(for recent high precision measurements of the transition point see 0). It is this order parameter, a monopole 
creation operator, that allows to rigorously distinguish between the Confined and the Coulomb phases of U(l)lgt. 
Other equivalent order parameters have been constructed later Q. 

In (1^ , a preliminary study of the mass of the monopole had allowed to claim that the DS exhibited in the Confined 
phase of U(l)lgt was of type II. In ^ similar study has been attempted with the Villain action but the authors 
could not gain enough statistics to study the monopole mass in the confined phase. 

After the work of Seiberg and Witten [l8j supporting the DS picture of confinement in the context of supersymmetric 
N = 1 gauge theories, we proposed an effective field theory (EFT) describing low energy physics for 4D U(l)lgt in 
terms of monopole fields [l2|. The EFT should work regardless of the action chosen on the lattice provided the 
correlation length is sufficiently large. 

In a later work, we tried to compare the predictions of the spectrum extracted from the EFT with the ones available 
from lattice simulations. This revealed some inconsistencies in lattice spectra [l^. Relying on a theorem in [ij we 
tried to match the results of with the ones of (22L 123^ (in the latter case the comparison can only be qualitative 
as these results are obtained with Villain action). This theorem states that, in the confined phase of 4D U(l)lgt, 
the Hilbert space of magnetically charged particles is contained into the Hilbert space of the neutral ones. Masses 
extracted from operators creating states with the same quantum numbers (angular momentum, parity and charge 
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conjugation (J^^)) but with different magnetic charges should coincide. The monopole mass of should coincide 
with the 0"*""*" gaugeball mass of and the dual photon mass of fl6| should coincide with the 1 gaugeball mass 

of 113. 

This, as we already pointed out was clearly not the case even at a qualitative level. In [l^l, the monopole was 
heavier than the dual photon whereas in [2^l23| the scalar gaugeball was lighter than the axial vector gaugeball. 

Later on, new results were obtained for the dual photon mass studying the electric flux tube profile between static 
charges |2J, |25j . The results presented therein (assuming the dual superconductor picture) reconcile the dual photon 
mass with the one of the axial vector gaugeball. The only piece of the spectrum in disagreement with the theorem 
was the monopole mass as calculated in jllf| . 

Our main result is to reconcile the value of the mass extracted from monopole correlators with the one of the 0^^ 
gaugeball. 

For that purpose, we introduce a modified order parameter as explained in the next section. Then, we present 
the detailed analysis of data obtained simulating this new order parameter on the lattice. We check that a single 
particle is responsible for the decay of its two-point correlation function. We extract the mass of this particle and find 
complete agreement with that of the 0++ gaugeball. Our method allows to extract a very good signal to noise ratio 
for the masses with an amount of data one order of magnitude smaller than that needed for existing techniques. We 
conclude with a discussion of the physical scenario emerging from our study. 

THE OPERATOR FOR SPECTRAL STUDIES 

The partition function of 4D U(l)lgt is defined as: 



z{(3) = / ( n n ) (-^^) ■ 



, (n,t) M=0 



We chose to consider the Wilson action: 



S = -Y. icosd0,^in,t)-l). (1) 

{n,t) i>fj.=0 

where {n, t) are the link variables of the four dimensional lattice whose sites are labelled by (n, t) . We use i as 
index for spatial directions and fi as index that spans all the four directions of the lattice. The field 9fj_{n,t) take 
value in C/(l) and Y[(n t) nji=o -^^mC""' the Lesbegue measure for each variable. We will abbreviate it as (H^^)- 
dOif^{n, t) is the lattice field strength term obtained acting with the exterior derivative on the link variables. d9i^{n, t) 
is hence defined on plaquettes identified by the coordinates {n,t) lying on the plane i — ji. At last, (3 — ^ with g 
being the U{1) the coupling constant. 

The order parameter for the Coulomb-Confined transition we consider was introduced by the Pisa group 
is the mean value of a monopole creation operator and shifts the plaquette field strength at a given Euclidean time by 
the contribution of the vector potential produced by a static magnetic source b{x). In the continuum, using standard 
notation this would be: 

fi{y,t)\A{x,t)) = \A{x,t) + -b{x-y)} (2) 

e 



with 

lj{y, t) = exp 

Further details on the subject can be found directly in [l^ . 
The lattice version of the operator is defined as: 



i- / d?xE{x,t)b{x-y) 



(3) 



^i{y,t) = exp^^/3(cos(6^(2/ - n) - d0io{n,t)) - cos{d9M{n,t))) (4) 
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If we consider the effect of a monopole antimonopole pair placed at points {x,ti) and (y, ^2) on the lattice the 
corresponding correlator will be : 




3 

Yl J|exp/3(cosd0,o(fi,O"l) (5) 

n,t<f{ti,t2}i=i 
3 

In this formula hi{x, y, n) — hi{x — n)t=ti — bi{y — n)t=t2 is the i component of the lattice vector potential produced by 
the monopole anti- monopole pair. It has support only on the two time slices ti and ^2- From now on we will consider 
the case for which x — y and the external field will be denoted bi{x — n). 

The expectation value (@J) is very difficult to extract from Monte-Carlo simulations. In order to see why this is so 
we define a new operator 



0{At)= Yl Y[e^ppicos{b,{x-n)-d0M{n,t))-l) (6) 

n,te{ti,t2} i=l 

where, to simplify the notation, we have ignored its dependency on all variables but At = 1^2 ~ ^i| • This is justified 
by the fact that we will always deal with expectation values of © that really depends (among other variables) on At. 
We introduce a modified action: 

3 3 
S = -Y. {cosdd,^{n,t)-l)- 

It is important to notice that the action iQ differs from the standard action ^ only on the two time-slices ti and t2 
where the monopole operator has support. On these time-slices Q vanishes. 

We can express the two point correlation function (jSJ as the mean value of the operator © on configurations 
generated with Q: 

{fi{x, t,)fi{x, t2)) = I / (n {-(^^) ^(^^) 

In this way the difficulty in measuring ^ appears immediately: we are trying to extract the mean value of © on 
configurations sampled uniformly on its support. On the time-slice ti and t2, support of JB)), the term — cos{d9io{n, t)) 
coming from Q cancels the analog contribution coming from the action (Q. In this way, on these time slices, 
importance sampling is lost. 

This is one of the reasons that forced the authors of 0| to introduce: 

d 

P = --g^^og{ii{x,ti)^i{y,t2)) (9) 
With this definition one recovers the importance sampling on the time-slices ti and ^2- On them the action becomes: 

3 

S'{t,,t2)= ^cos(&,(f-n)-d0,o(ri,t))-l (10) 

n,t(E{ti,t2} i=l 

If one is interested in studying the mean value of p this is enough In our case we need to go a little further as 
we want to extract masses in the confined phase. To do this we are forced to study the decay of (O as function of 
the time separation t between monopole and antimonople. Assuming that for large t the decay of (O is driven by a 
single particle, in the confined phase, we obtain: 

{^l{x, t)fl{x, 0)) AAf i/2t-3/2e-^^*. (II) 
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In this formula /z is the v.e.v of the monopole operator ij,{x,t). It is different from zero in the confined phase. A is 
the projection ofthe monopole state on the vacuum . M is the monopole mass. In order to keep formulas compact 
we are neglecting the effect of periodic boundary conditions but we will add them when studying numerical results. 

Fluctuations of in the confined phase, are caused by fluctuation of the v.e.v. /i^ and completely screen the 
decay we need to unmask to extract M. To deal with this in a very huge amount of data was collected. We 
can, however, introduce a new definition of p that, following the ideas introduced in in the context of gaugeball 
spectroscopy, approximates the connected part of (jsj and eliminates p?. The new operator is defined as: 

p' = ^log(/x(f,t)A2(f,0)). (12) 

With this definition importance sampling on the monopole anti-monopole time slices is kept (as in the case of 
and fluctuations at scale are eliminated. Starting from Hll|l we obtain the large t behaviour of H12|l : 

p' ^ ^Mi/2t-3/2g-A« _ (^^^2 ^ AM'/h-'/^e'^") ~' (13) 

This expression has 3 free parameters: /i, M, A. On a large enough lattice we can expand H13() in a regime where 

2tJ ^2 

from which it is easier to extract a precise determination of the mass. However the results presented in this paper are 
obtained with the full expression H13() . 

The simulation algorithm for studying p' decays is very similar to the one needed to study p decays. Starting from 
the expression (jH)) and taking the logarithm we get 

dt log < p{x, t)fi{x, 0) ^ ^ (14) 

We need to use the first order approximation in dtf'. 

dtef^'^\t^to=ef^'Hdtf)\t^t, (15) 
In this way we can express the derivative as : 

3 

dtO{t) = 0{t) ^ 51 ^* (/3(cos(6,(n) - de,o{fi, to)) - 1) 

n i=l 

and we can include the factor 0{t) in the measure. This defines a new action: 

3 

S' = S + ^Y^ cos(6,(n) - d9M{n, t)) - 1. 

n i=l 

The expression p4f) in this way becomes: 

dt\og{p{x,t)p{x,0)) = (J[ve) e-P'''dt{logO{t)-(3S{t)). (16) 

where S{t) is the Wilson action defined in Q restricted on the time shces t and 0. We hence need to compute the 
mean value of 

dt{\ogO{t)-ps{t)) 

on configurations generated with S': 

p' ^{dt{logO{t)-pSm' (17) 
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FIG. 1: This figure represents the configurations we used to calculate the discrete derivative. The circle represents the time 
direction of our lattice. The thicks represent monopole and anti-monopole terms included in the action. The solid points 
represent the time slices used to measure the operator and construct the derivatives. On the left hand side we see the 
symmetries of the construction when interchanging monopole and antimonopole derivatives. On the right hand side differences 
between symmetric and backward derivatives are sketched. 



To calculate the derivative of logC'(t) — PS{t) numerically in principle we could choose the forward, backward or 
symmetric prescriptions. We are considering, however, the approximation H15|) with f{t) = dt (logC'(i) — S{t)) and 
hence need to minimize f(t). This requirement immediately selects the symmetric prescription f'(t) — {f{t + 1) — 
f{t — l))/2. The reason for that is outlined in the right part of figure There, we sketch the scenario used in 
our computations. The circle represents the temporal direction of the lattice (the other directions are omitted for 
simplicity) and the dashed and solid thicks represent the insertion of monopole and anti-monopole operators in the 
action (|10|l at the two time slices separated by a distance t. The solid dots represent the time slices used to calculate 
the derivative in the symmetric case and in the back-ward case. Symmetric derivatives involve time slices separated by 
two lattice spacings. Nevertheless, their plaquettes are updated with the same action and are hence more correlated 
than the adjacent time slices involved in the backward or forward derivatives calculations. These last are updated 
with different actions since on one of them there is the monopole (or antimonopole) contribution. The symmetric 
prescription, hence, minimizes dt (logO(t) — S{t)) and is used to extract the correct value of p' . 

In order to obtain H12|l from (|SJ) we have a further freedom: we can derive the correlator either with respect to the 
monopole or to the antimonopole position, the only difference being an overall minus sign. This is shown on the left 
part of figure n Derivatives on the monopole field amount to minus derivatives on the anti-monopole field. 



NUMERICAL RESULTS 

We performed three runs of simulations on three different lattices. The smallest is a 6^ x 12 lattice the medium 
is a 8^ X 16 and the largest 10^ x 20. We collected over 400000 sweeps of statistics for all the lattices . We used the 
standard 1 heat-bath and 3 over-relaxation algorithm. In order to extract the value of the masses from the numerical 
data, we have to include the effect of periodic boundary conditions (p.b.c). Once we define f{x) = AM^/^x^'^^^e'^^^ , 
the p.b.c. modify Hll|l producing: 

{^^{x, ii)/i(f , t2)) ^ A*' + f{r) + fin) + 3 (/(r,) + /(r,*) + /(r,,) + f{rsst)) ■ (18) 

where rtrs,rss take into account the different path we can chose to go from {x,ti) to (x, t2). We can in fact chose 
the shortest path (r = t = \t2 — ti\) or the path that winds once in the temporal direction (r^ — L — t); we can 
also wind once around the spatial direction (r^ ~ ^/t'^ + L"^, we have three different choices), twice around the 
spatial direction (r^s = ^yW+{2Ls)^ we have three possible choices) once around the spatial and the time direction 
(rst = ^ (L ~ t)^ + Li we have three possible choices) and twice around the spatial direction and one around the time 
direction (rsst = \/ [L — ty + (2Ls)^ we have three possible choices). Taking into account all these possibilities we 
extracted the masses depicted in figureOlin a range of /3 which ensured the system is in the confined phase. The first 
thing we notice form figure El is that, for all the values of f3 considered, the /n.d.f. obtained fitting the numerical 
data with expression l|13|l is lower than one. The only exceptions are, indeed, points very close to the transition 
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/n.d.f . 

10 X 10x10x20 



O Anlimunopule (lerivalive 
□ From monopole derivalivt 



FIG. 2: This figure is a plot of the /n.d.f . for the fits we performed with expression l|18^ on the lattice 10'' x 20. For each /3 
we have at least one channel with /n.d.f . <g 1. This confirms the validity of (1131 to describe the correlator decay. 



(which is first order with very long correlation length) and are surely due to metastabilities caused by our updating 
algorithm. This means that the leading contribution to the decay is correctly described by a single particle excitation. 
As explained in the previous section, we use the symmetric prescription for the derivative in expression H17|l and 
derive both with respect to the monopole and anti- monopole position . From the figure |21 it is clear that the masses 
extracted in both cases are compatible within error-bars. 

In figure 0] we show typical data-points we obtain for the derivative of the correlation function (|17|l from the 
simulations performed. The line represents the best-fit curve. The upper plot contains points obtained deriving 
with respect to the monopole position whether the lower plot contains points obtained deriving with respect to the 
antimonopole position. In both cases we where forced to add a constant term to the expression (|13l) to correct from 
systematic errors induced by discretization and by expanding the exponential of the derivative in H15|l . The best-fit 
curves are obtained using a subset of distances (from 4 to 16) to avoid contamination from higher states in the same 
channel and safely consider only the single particle decay. 

Comparing our results with the ones in |l6l | one notices a clear improvement of the signal to noise ratio. The 
difference between the values of the masses we extract and the ones obtained in 0| is probably due to the high 
noise that prevented the authors of from using the full expression H18() in fits. The authors, forced to use a zero 
momentum approximation of the decay Hll|l . introduced a systematic error that shifted the masses from their actual 
value. 

We can also safely identify the particle responsible for the decay. The masses we get from the fit with the complete 
expression H18f) are indeed fully compatible with the known value for the gaugeball O"'""'' mass (see i.e. for a recent 
high precision study). To show this we took into account finite size effects in the two plots |S| 

In the first we considered the results we obtained for the masses (from both monopole and anti-monopole channels) at 

(3 = 1.005. We plotted them versus the lattice size (calculated as L = [L^ x Lt)^^*). The last point at L = 16 is taken 
from the results of the 0+"'" mass contained in Under the assumption that the discrepancy between our results 
and the one |^ is due to finite size effects we made a two parameters fit with the expression M(L) = Af^o — AL~^. 
In the second plot there is the same study at (3 = 1.0099. The values safely confirm the identification of the 
particle responsible for the decay Hill) with the 0++ gaugeball. Furthermore the fact that the thermodynamic limit 
is approached from below is the expected behaviour in the confined phase (see i.e. H^). These results confirm the 
validity of a theorem about the structure of the Hilbert space in the confined phase of compact lattice U{1) in four 
dimensions proved in there are no super-selected sectors labelled by the magnetic charge. The last point to 

stress is that our method prevents a precise measurement of the v.e.v. (as it is designed to get rid of it as far as 
possible). This is shown in figure El and implies that, if interested in studying the order parameter for the Coulomb 
Confined transition, one should better use the standard p definition (isf. 
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Monopole masses in the confined pliase 

8x8 X 8 X 16 lattice and 400000 sweeps 
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Monopole masses in the confined phase 

10 X 10 X 10 X 20 lattice and 400000 sweeps 



Simm. der. on mon. pos. 
^ simm. der. on Amon. pos. 



0.8 -■ 

0.6 - 
0.4 - 



f i I 



0.2 



FIG. 3: These figures collect masses calculated with the two possible strategies on the three lattice sizes we considered in a 
range of /3 in the confined phase. We can derive the correlation function either with respect to the monopole position or to the 
anti-monopole position. Masses extracted with any of these prescriptions are compatible within error bars. 



CONCLUSIONS 



In this work we considered the spectrum extracted from magnetically charged operators. We introduced a new 
powerful technique that reduces, with respect to existing methods, one order of magnitude the amount of data 
necessary to extract a good signal for correlations of such operators . 

The improvement of the technique relies on the fact that it allows to study the connected part of the correlation 
functions. As a first application we study the Hilbert space of confined compact four dimensional C/(l) theory. We 
consider monopole anti-monopole correlation functions on three lattice sizes (6'^ x 12, 8'^ x 16, 10'^ x 20). We check 
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FIG. 4: This figure represents the correlation function decay for a 10^ x 20 lattice &t {3 — 1. The points on the upper plot 
are obtained deriving with respect to the monopole field, while the point on the lower part deriving with respect to the 
anti-monopole field . 



that the large time correlation functions decay is driven by a single particle. We also identify this particle by studying 
its mass. The value for the mass we extract is compatible with the one found in literature for the gaugeball 0"'"'". This 
is in complete agreement with a theorem proved in jij stating that the Hilbert space of compact lattice U{1) in four 
dimensions does not contains any magnetically super-selected sectors. Completing our results with the ones obtained 
by other groups with different techniques |22,|2J,|25j one can have a full picture of the Hilbert space of compact U{1) 
lattice gauge theory in four dimensions as the space spanned by gaugeball states. We are working on similar studies 
for the case of non Abelian lattice gauge theories. We are also working on the implementation of an new algorithm 
to measure directly the order parameter t) and its finite size scaling based on ideas similar to the ones presented 
in this work. 
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